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I. 


INTRODUCTION 


Law  Of  Revelation :  "  The  Hidden  Flew  Never  Remeins  Hidden .  " 


Simply  stated,  the  dynamically  aimed  free-f light  rocket  (DAFR  or  DAFFR) 
concept  is  to  decrease  the  sensitivity  of  an  artillery  rocket  to  other  error 
sources  by  increasing  wind  sensitivity.  The  wind  effect  is  reduced  during 
firing  by  round-to-round  correction  in  aiming  based  upon  measurements  of  the 
previously  fired  rounds.  The  concept  has  been  described  in  detail  by  McCorkle 
and  Lilly  1 ,  and  some  refinements  to  the  error  analysis  were  later  made  by 
Gibson.2 

It  was  assumed  1»2  that  the  measurements,^  z  ,  would  be  averaged  to 
obtain  an  unbiased  minimum  variance  estimate,  x  , 


A 

X 


(1) 


The  variance  of  the  estimate  is 


2 

°x 


(2) 


2 

and,  assuming  the  variance, 


is  constant. 


2 


(3) 


The  variance  of  the  wind  is  not  constant  and  was  modeled  as 


a 


2 

u 


N2  yz 


(4) 


where  N  is  the  number  of  rounds  fired  and  y  is  the  sensitivity  slope  which 


1  William  C.  McCorkle,  Jr.,  and  J.  A.  Lilly,  "An  Adjusted  Fire  Technique 
for  a  Highly  Accurate  Free  Flight  Rocket  Artillery  System,"  US  Army  Missile 
Command,  Tech.  Rpt.  RD-74-13,  Redstone  Arsenal,  AL,  25  June  1974  (AD  B007539L). 

2  J.  D.  Gibson,  "Error  Analysis  of  a  Dynamically  Aimed  Free  Rocket," 
Texas  A&M  University,  College  Station,  TX  (Unnumbered  and  undated,  not  held 
by  OTIC,  formerly  DDC,  or  RSIC). 


was  derived  from  data  in  the  Tactical  Artillery  Rocket  Environment  (TARE) 
report. 3 

Since 


N(N  +  1)  (2N  +  1) 
6 


the  resulting  variance  is 


(5) 


2 

ax 


(l  +  I\  02  +  + 

'  .  (N  +  1)(2N  +  1)1 

V  nJ  ap  N 

.  6N  J 

(6) 


where  a2  is  the  sum  of  the  squares  of  all  Gaussian  white  plant  noises,  and 
a*  is  the  variance  of  the  observation  noise.  Equation  (6)  eventually  diver¬ 
ges  as  N  becomes  large  because  the  linearity  assumption.  Equation  (4),  breaks 
down.  Though  additional  data  may  not  improve  an  estimate,  it  should  not  make 
the  estimate  worse. 

The  TARE  report3  states  that  "Readings  were  taken  with  aerovane  type 
anemometers  located  at  heights  of  approximately  60  feet  and  using  a  smoothing 
circuit  of  one  (1)  minute."  The  times  between  rounds  previously  considered1 »2 
were  around  five  to  ten  seconds.  The  time  delay  characteristics  previously 
used  were  most  likely  those  of  the  "1  minute  time  filter"  not  the  wind! 

This  leaves  us  with  two  questions: 

•  What  are  the  (unfiltered)  characteristics  of  the  wind? 

•  What  is  their  effect  on  the  error  analysis? 


II.  WIND  AND  TURBULENCE 


Thumb 1  S  First  Postul  ate :  "It  is  better  to  solve  a  problem  with  a  crude  approximation 
and  know  the  truth,  f  10%  ,  than  to  demand  an  exact  solution  and  not  know  the  truth  at  all!  " 

The  wind  model  presented  below  was  developed  mainly  from  the  results  of 
a  literature  search  on  atmospheric  turbulence  done  by  Stewart.4 


3  Tactical  Artillery  Rocket  Environment  (TARE)  Committee,  "Tactical 
Environment  for  Large  Free  Flight  Rocket  Systems,"  US  Army  Ballistic  Missile 
Agency  Report,  R-S-61-1,  Redstone  Arsenal,  Alabama,  1  June  1961 

(AD  B952069L) . 

4  D.  A.  Stewart,  "A  Survey  of  Atmospheric  Turbulence  Characteristics," 
US  Army  Missile  Command,  Technical  Report  RR-81-6,  Redstone  Arsenal,  Alabama, 
19  August  1981. 
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Figure  1  is  a  plot  of  the  horizontal  wind  power  spectral  density.5*6  The 
spectral  gap  between  eight  hours  and  ten  minutes  serves  to  explain  why  we  may 
speak  of  an  average  wind.  To  reduce  randomness  one  should  operate  as  far  as 
possible  on  the  other  side  of  the  peak  due  to  turbulence. 

Assuming7’4  that  the  longitudinal  wind,  u  ,  is  the  sum  of  the  longitudi¬ 
nal  average  wind,  u  ,  and  the  turbulent  wind,  u'  ,  then 

u ( t )  =  u  +  u "(t)  (7) 


by 


7  4 The  turbulent  wind  is  related  to  the  turbulent  wind  at  some  previous  time 


u'(t  +  t)  =  p(r)  u"(t)  +  u " ( t )  (8) 

whr"e  p(t)  is  the  autocorrelation  coefficient  for  a  time  delay,  t  ,  and 
u"U)  is  the  random  component  of  the  turbulence. 

The  variance,  a7,,  ,  of  the  random  components  is  defined  by  the  relation¬ 
ship7’4 


au"  =  (9) 

so  that  the  turbulent  energy,  a7.  ,  be  conserved  with  time. 

The  standard  deviation  due  to  turbulence,  ,  is  related  to  the  average 
wind,  U"  ,  by  the  approximation6’4 


a 


u ' 


(10) 


where  h  is  the  height  and  zQ  is  the  roughness  length8.  Figure  2.  The 


5  A.  J.  McDonald,  Wind  Loading  on  Buildings,  New  York:  John  Wiley  &  Sons, 
1975. 

6  J.  L.  Lumley  and  H.  A.  Panofsky,  The  Structure  of  Atmospheric  Turbulence, 
New  York:  John  Wiley  &  Sons,  1964. 

7  S.  R.  Hanna,  "Some  Statistics  of  Lagnangian  and  Eulerian  Wind  Fluctu¬ 
ations,"  Journal  of  Applied  Meteorology,  Vol  18,  April  1979,  pp.  518-525. 

8  J.  W.  Kaufman  (ed.),  "Terrestrial  Environment  (Climate)  Criteria  Guide¬ 
lines  for  Use  in  Aerospace  Vehicle  Development,"  1977  Revision,  Second  Edition, 
NASA  Technical  Memorandum  78118,  June  1979. 


Figure  1.  Power  spectrum  of  horizontal  wind  speed. 
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Figure  2.  Roughness  length  for  various  terrains. 


lateral  components  are  approximately9 


av  =  (0.64)ou. 


(ID 


and 


s,  -  ‘°-52)v 


(12) 


Of  course,  "v  =  w  =  0  by  definition.  These  approximations  hold  best  for  a 
neutral  atmosphere,  which  implies  strong  winds. 

Figure  3  is  a  plot  of  the  intensity,  a  ^/t r  ,  against  height  for  various 
roughness  lengths. 

The  Dryden  correlation  coefficients9  are  as  follow: 

•  Longitudinal 


PU(T) 


-u  x/L, 


(13) 


•  Transverse 


Pw(t)  =  e  V(1  -  ux/2L  ) 


pw(t)  =  pv(t) 


(14) 

(15) 


where  L  is  the  isotropic  turbulence  integral  scale. 
In  the  atmospheric  boundary  layer9. 


25  h 

,0.4 


(16) 


where 


c  =  e 


-0.025(£n  zq)  +  0.17  £n  zQ  -  0.8 


(17) 


9  W.  Frost,  B.  H.  Long,  and  R.  E.  Turner,  "Engineering  Handbook  on  the 
Atmospheric  Environmental  Gudielines  for  Use  in  Wind  Turbine  Generator  Devel¬ 
opment,"  NASA  Technical  Paper  1359,  December  1978. 


8 


and 


Ly  =  =  (0.4)h  ,  h  <  250  meters.  (18) 


For  h  >  250  meters 


Lu  =  Ly  =  Lw  =  250  meters  (19) 


Figure  4  is  a  plot  of  l_u  against  height  for  various  roughness  lengths.  To 

illustrate  how  the  formula/figures  may  be  used,  consider  a  height  of  5  meters 
and  a  roughness  length  of  0.1  meter.  From  Figure  3 

ou.  =  (0.25)  u  (20) 


and  from  Figure  4 

Lu  =  100  meters  (21) 


Substituting  into  Equation  (13)  and  then  Equation  (9)  results  in 

au„  =  (0.25)  u  (l  -  e'2*17100)  (22) 

A  plot  of  Equation  (22)  is  shown  in  Figure  5.  L  and  L  are  so  short, 

V  W 

Equation  (18),  that  the  transverse  components  are  essentially  white  noise  and 
correction  is  not  feasible,  though  they  would  contribute  to  the  error  budget. 

The  above  model  for  the  turbulent  wind  corresponds  to  a  Gauss-Markov 
process,  that  is,  Gaussian  white  noise  passed  through  a  low  pass  filter  where 
the  order  of  the  process  is  the  order  of  the  filter10.  This  model  is  only 
reasonable  for  the  most  right-hand  portion  of  Figure  1,  but  for  time  delays  of 
less  than  one  minute  it  is  within  20  percent  of  observed  autocorrelation 
coefficients.7  An  in-depth  discussion  of  the  bases  for  the  model  will  be 
found  in  the  report  by  Stewart.4 


l0A.  Gelb  (ed.).  Applied  Optimal  Estimation,  Cambridge:  MIT  Press,  1974. 
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Figure  4.  Turbulence  integral  scale  versus  height. 


c/uiOZ 


Figure  5.  Turbulence  standard  deviation  versus  time  delay. 


III. 


THE  ERROR  ANALYSIS  -  REVISITED 


Grossman's  Misquote  of  H.  L.  Menchen: 

"Complex  problems  hove  simple,  eesy  to  understood  wrong  answers." 

Instead  of  using 

*  N^2 

use  the  following  in  Equation  (2): 


,  2  (, 

«;■  ■  oj-  u  -  •  / 


Since 


k=0 


-ak  _  1  -ea(N  +  11 

'  '  “IT 


for 


x  =  N  T 


2 


w  E  (v  *  +  °u“2) 


k=l 


(' 


(e: 

(2: 

(2< 

(2! 


becomes 


2  2  2 
-  ,  _  °P  +  %  .  V 

°e  '  N - +  ~W 


N  +  1 


1-e 


■2(N+1)IT  T/L 


.1-e 


-2u  T/L, 


(21 


For  large  u  T/Lu  ,  Equation  (26)  approaches 


n  2  +  a  2  +  a  2 


'  2 

'  N 


P  V  u 


(2; 


that  is,  the  Gauss-Markov  process  "becomes"  essentially  Gaussian  white  noise 


Since 


i. 


2 


+  a., 


+ 


(28) 


substituting  Equation  (26)  yields 


+  au 


2 


-2(N+l)uT/L  \ 
l-e  u  \ 

- ^uT7l - ] 

>  l-e  u  /J 


(29) 


Plots  of  Equation  (29)  are  given  in  Figures  6  through  11.  The  standard 
deviations  have  been  arbitrarily  expressed  in  mils;  though  representative,  the 
values  do  not  reflect  an  actual  system.  The  quantity,  uT/l_u  »  requires  some 

explanation.  For  a  given  turbulence  integral  scale,  for  example. 


Lu  *  100  meters 


and  an  average  wind 

IT  =  5  meters  per  second 


for 


u 


it  follows  that  the  time  between  rounds  must  be 
T  =  2  seconds 

For  a  small  precision  error.  Figure  6,  the  improvements  from  decreasing 
the  time  between  rounds,  T  ,  are  dramatic.  For  a  large  precision  error. 
Figure  11,  there  is  not  much  improvement  at  all  since  that  error  dominates. 

These  figures  may  be  used  for  other  integral  scales,  Lu  ,  and  average 
winds,  u  ,  but  the  time  between  rounds,  T  ,  would  need  to  be  recomputed. 
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Figure  6.  Number  of  round*  already  fired. 


Figure  7.  Number  of  rounds  already  fired. 


Figure  8.  Number  of  rounds  already  fired. 


Figure  9.  Number  of  rounds  already  fired. 


10  mH 


Figure  10.  Number  of  rounds  already  fired. 
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Figure  11.  Number  of  rounds  already  fired. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Wetheru's  Law  of  Suspended  Judgement: 


"Assumption  is  the  mother  of  all  screw-ups," 


The  results  of  incorporation  of  the  revised  wind  (turbulence)  model  in 
the  error  analysis  supports  the  overall  conclusions  of  the  previous  studies.1’2 
Although  it  would  appear  that  the  time  between  rounds  should  be  shorter,  it 
does  place  those  conclusions  on  a  bit  firmer  foundation. 


Probably  the  most  significant  contribution  of  this  effort  has  been  the 
incorporation  of  surface  roughness  via  zq  .  Fichtl  and  McVehil11  report  a 

roughness  of  around  0.2  meter,  where  downwind  of  vegetation  1  to  2  meters  in 
height,  but  0.3  meter  when  450  meters  downwind,  and  0.4  to  0.6  meter  when 
200  meters  downwind  of  trees  10  to  15  meters  tall.  It  would  appear  desirable 
not  to  fire  from  positions  downwind  of  trees  or  buildings  if  the  tactical  sit¬ 
uation  permits.  Deaves12’13  has  proposed  a  computational  wind  model  which  may 
be  of  use  in  analyzing  the  effects  downwind  of  changes  in  surface  roughness. 


Of  course,  if  the  average  is  taken  over  too  long  a  time  one  would  have 

IT  ,  not  u ( t+T )  .  If  the  DAFR  concept  proceeds  into  development  after  demon¬ 
stration,  it  is  recommended  that  a  probabilistic  rather  than  a  statistical 
filter  be  developed,  that  is,  a  Kalman  filter10  be  utilized.  The  above  turbu¬ 
lence  model  is  suitable  for  Kalman  filter  formulation 


(30) 


for  the  state  equation,  and 


2k+l  " 


k+1 


(31) 


11  6.  H.  Fichtl  and  G.  E.  McVehil,  "Longitudinal  and  Lateral  Spectra  of 
Turbulence  in  the  Atmospheric  Boundary  at  the  Kennedy  Space  Center,"  Journal  of 
Applied  Meteorology,  Vol  9,  February  1970,  pp.  51-63. 


12  D.  M.  Deaves,  "Computations  of  Wind  Flow  Over  Changes  in  Surface  Rough¬ 
ness,"  Journal  of  Wind  Engineering  and  Industrial  Aerodynamics,  Vol  7,  1981, 
pp.  65-94. 


13  0.  M.  Deaves,  "Terrain  -  Dependence  of  Longitudinal  R.M.S.  Velocities 
in  the  Neutral  Atmosphere,"  Journal  of  Wind  Engineering  and  Industrial  Aero¬ 
dynamics,  Vol  8,  1981,  pp.  259-275. 


i 


[  for  the  observation  equation,  where  is  just  u^  with  observation 

|  noise,  ,  added.  Of  course,  these  state  and  observation  equations  would 

[  have  to  be  expanded  to  include  the  missile  and  tracker  dynamics.  Besides 

[  explicit  incorporation  of  the  Markov  process  in  the  filter,  an  error  analysis 

I  would  be  generated  in  the  filtering  process. 

L 

l 

fe  There  is  no  lew  except  the  law  that  there  is  no  law." 


-John  Archibald  Wheeler 
Physicist 
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